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Abstract. 



The semileptonic decay of the neutral K meson, K\ — > 7r ± e =F z/ (K e3 ), was used to study 
the strangeness-changing weak interaction of hadrons. A sample of 5.6 million reconstructed 
events recorded by the NA48 experiment was used to measure the Dalitz plot density. Ad- 
mitting all possible Lorentz-covariant couplings, the form factors for vector (f + (q 2 )), scalar 
(fs) and tensor (/r) interactions were measured. The linear slope of the vector form factor 
A+ = 0.0284 ±0.0007 ±0.0013 and values for the ratios |/ 5 //+(0)| = 0.015lg;g?J ±0.012 and 
|/t//+(0)| = 0.05^g;o4 ± 0.03 were obtained. The values for f s and fx are consistent with 
zero. Assuming only Vector-Axial vector couplings, A + = 0.0288 ± 0.0004 ± 0.0011 and a 
good fit consistent with pure V-A couplings were obtained. Alternatively, a fit to a dipole 
form factor yields a pole mass of M = 859 ± 18 MeV, consistent with the K*(892) mass. 



1 Introduction 

The study of semileptonic decays of K L mesons gives valuable information about the 
strangeness-changing weak interaction and can be used as a test for possible non-vectorial 
components of the weak hadronic interaction and of models of low-energy strong interactions. 

The most general form of the matrix element for K e3 decays is given by[l] : 

M = sJl/2G F sin e c u- u (l - 75 ) [m K f s + (1) 
(l/2i){{P K ± P n )J + ± (P K - P«) x f-}lx + i(fT/m K )a XT [P K ] x [P v ] T ] Ul , 

where 9 C is the Cabibbo angle, and ui are the lepton spinors, Pk and P n are the kaon and 
pion four-momenta respectively and ttik is the kaon mass. The determination of the scalar 
(fs), vector (/ + and /_) and tensor (fa) form factors is based on measurement of the Dalitz 
plot density which in the kaon rest frame has the following form [2]: 

P (K, K) - *lh(q 2 )] 2 + c[f s + —(e;- Ei)f T }\ (2) 

rriK 

where E* is the energy of the particle i in the kaon rest frame, 
f + (q 2 ) = /+(0)(1 ± \ + q 2 /ml + A' + g 4 /m 4 J, 
a = m K (2E* e E;-m K E 7T f ), 
c = m 2 K E n ', 

T? I _ ( m K+ m l) _ 77* 

n 2m K 1"' 

q 2 = (m\ +ml- 2m K E* n ) 

In the expression for the vector form factor V, the /_ contribution can be neglected because 
it is proportional to the electron mass squared. The dependence of the / + form factor on 
the momentum transfer q 2 and its value /+(0) at q 2 = are of theoretical interest [3], [4], 
[5]. The q 2 dependence is usually assumed to follow a pole-dominance formula f+(q 2 ) = 
/+(0)/(l - q 2 /M 2 ), where M is the K* meson mass, M = 892 MeV. This leads in a linear 
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approximation to A + = rn^/M 2 = 0.0245 if one uses the 7r + mass for m n . A calculation 
in chiral perturbation theory to order 0(j» 6 ) [3] predicts a value of A + = 0.022. Other 
calculations in chiral perturbation theory [4] and in lattice QCD [5] concentrate on the value 
of /+(0), which is important for the absolute rate of K e3 decay. 

Evidence for non-zero scalar and tensor form factors in the case of K + — > 7i e + u has been 
reported in [6]. Recent measurements [7], [8] of the charged decay modes have not confirmed 
the results of [6] however, while investigations of neutral kaon decays have revealed no signif- 
icant deviation from vector type interactions ([9], [10]). The study presented here improves 
the statistical and systematic significance of these investigations using a measurement per- 
formed with the NA48 detector in a neutral kaon beam at the CERN SPS. 



2 Experimental setup 



The NA48 detector was designed for a measurement of direct CP violation in the K° system. 
Here we use data from a dedicated run in September 1999 where a Kl beam was produced 
by 450 GeV/c protons from the CERN SPS incident on a beryllium target. The decay region 
is located 120 m from the Kl target after three collimators and sweeping magnets. It is 
contained in an evacuated tube, 90 m long, terminated by a thin (3- 10~ 3 A ) kevlar window. 

The detector components relevant for this measurement include the following: 

The magnetic spectrometer is designed to measure the momentum of charged particles 
with high precision. The momentum resolution is given by 

= (0.48 ©0.009 -p)% (3) 

V 

where p is in GeV/c. The spectrometer consists of four drift chambers (DCH), each with 8 
planes of sense wires oriented along the projections x,u,y,v, each one rotated by 45 degrees 
with respect to the previous one. The spatial resolution achieved per projection is 100/xm and 
the time resolution is 0.7 ns. The volume between the chambers is filled with helium, near 
atmospheric pressure. The spectrometer magnet is a dipole with a field integral of 0.85 Tm 
and is placed after the first two chambers. The distance between the first and last chamber 
is 21.8 m. 

The hodoscope is placed downstream of the last drift chamber. It consists of two planes 
of scintillators segmented in horizontal and vertical strips and arranged in four quadrants. 
The signals are used for a fast coincidence of two charged particles in the trigger. The time 
resolution from the hodoscope is 200 ps per track. 

The electromagnetic calorimeter (LKr) is a quasi-homogeneous calorimeter based on 
liquid krypton, with tower read out. The 13248 read-out cells have cross sections of 2 x 
2 cm 2 . The electrodes extend from the front to the back of the detector in a small angle 
accordion geometry. The LKr calorimeter measures the energies of the e ± and 7 quanta by 
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gathering the ionization from their electromagnetic showers. The energy resolution is : 

a(E) ,3.2 9.0 , 

^ = (-^©-©0.42)% (4) 

where E is in GeV, and the time resolution for showers with energy between 3 GeV and 100 
GeV is 500 ps. 

The muon veto system (MUV) consists of three planes of scintillator counters, shielded 
by iron walls of 80 cm thickness. It is used to reduce the Kl — > 7i ± fi T h' background. 

A more detailed description of the NA48 setup can be found elsewhere [11]. 



3 Data processing and Monte Carlo simulation 



3.1 Trigger and data taking 



The trigger was built in two levels. In the first level trigger, the presence of at least two 
hits in the hodoscope was required. In the second level trigger, the drift chamber hits were 
used to reconstruct tracks and at least one pair of tracks was required to have a transverse 
separation below 5 cm within 4.5 K$ decay lengths from the end of the last collimator. A 
first level control trigger requiring at least one hit in the hodoscope in coincidence with at 
least two tracks segments in the first DCH, downscaled 20 times, was also implemented for 
trigger efficiency calculations. 



3.2 Event selection 



The data sample consisted of about 2 TB of data from 100 million triggers, with approxi- 
mately equal amounts recorded with alternating spectrometer magnet polarities. The data 
were reconstructed and subjected to off-line filtering. The following selection criteria were 
applied to the reconstructed data to identify K e3 decays and to reject background, keeping 
in mind the main backgrounds to K e3 , which are Kl — > 7i ± fi T u (K^) and Kl — > n + n^n 

(K 3n y. 

- Each event was required to contain exactly two tracks, of opposite charge, and a recon- 
struced vertex in the decay region. To form a vertex, the closest distance of approach between 
these tracks had to be less than 3 cm. The decay region was defined by requirements that 
the vertex had to be between 6 and 34 m from the end of the last collimator and that the 
transverse distance between the vertex and the beam axis had to be less than 2 cm. These 
cuts were passed by 35 million events. 

- The time difference between the tracks was required to be less than 6 ns. To reject muons, 
only events with both tracks inside the detector acceptance and without in-time hits in 
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the MUV system were used. For the same reason only particles with a momentum larger 
than 10 GeV (p m i n ) were accepted. In order to allow a clear separation of pion and electron 
showers, we required the distance between the entry points of the two tracks at the front 
face of the LKr Calorimeter (Di kr ) to be larger than 25 cm. As a result 14 million events 
remained. 

- For the identification of electrons and pions, we use the ratio of the measured cluster energy, 
E, in the LKr calorimeter associated to a track to the momentum, p, of this track as measured 
in the magnetic spectrometer. The ratio E/p is shown for all tracks of the 14 million events 
in fig.l. For the selection of K e3 events, we require one track to have 0.93 < E/p < 1.10 
(electron) and the other track to have E/p < 0.90 (pion). 11.7 million events were accepted. 

tracks/0.01 (binsize) 

4 : 

3.5 : 

3 : 

2.5 : 

2 : 

1.5 : 



1 - 




E/p 

Fig. 1. Distribution of the ratio of the shower energy E reconstructed by the LKr and the momentum 
p reconstructed by the spectrometer. 

- In order to reduce background from K 3n decays, we required the quantity 

p ,2 _ (m 2 K - m\_ - mlo) 2 - 4(m^m^ + m\p\) 

4(pl + ml_) [b) 

to be less than — 0.004(GeV/c) 2 . In the equation above, p± is the transverse momentum of 
the two track system (assumed to consist of two charged pions) relative to the K\ flight 
direction and m + _ is the invariant mass of the charged system. The variable Pq 2 peaks at 
zero if the charged particles are pions from the decay K 37T . The cut removes (98.94 ± 0.03)% 
of K 3w decays and (1.03±0.02)% of K e3 decays as estimated with the Monte Carlo simulation 
(Sect. 3.4). After this cut, we were left with 11.4 million K e3 candidate events. 

The neutrino momentum in K e3 decays is not known and the kinematic reconstruction of 
the kaon momentum from the measured track momenta leads to a two-fold ambiguity in 
the reconstructed kaon momentum. In order to measure the kaon momentum spectrum, we 
selected events in which both solutions for the kaon momentum lie in the same bin of width 
8 GeV. These 4 • 10 5 events we call "diagonal events". 
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The last selection criteria were the following two requirements: each solution for the kaon 
energy had to be in the energy range (60,180) GeV and the reconstructed electron and pion 
momenta in the center of mass of the kaon had to lie in the kinematically allowed region of 
the Dalitz plot for K e3 decays. 

As a result of this selection, 5.6 • 10 6 fully reconstructed K e3 events were selected from the 
total sample. 



3.3 Background 

The main background to K e3 events arises from and K 3n decays. The background from 
K^s was reduced by removing events with an in-time muon signal and by the E/p require- 
ments. The probability for a pion from to fake an electron (E/p > 0.93) was measured 
to be 5.7- 10~ 3 . This was obtained from a sample of 75000 pion tracks from K e3 events where 
the other track had E/p > 1.02. The inefficiency of the MUV system for the total exposure 
was less than 2%. From this we estimated the number of fake K e3 events from K^ 3 to be at 
most 400 ± 100 events. For the suppression of background from K 3n decay, we used the cut 
Pq 2 < — 0.004(GeV/c) 2 and the E/p requirements. The number of background events from 
this source is estimated to be less than 20 events. 



3.4 Monte Carlo Simulation 

The NA48 detector response was simulated using a GEANT-based Monte Carlo (MC) [11]. 
The kaon energy spectrum in the MC was taken from the kaon energy distribution recon- 
structed using diagonal K e3 events, as defined above (Sect. 3.2). 

For the generation of K e3 MC events, the slope parameter A + and the scalar and tensor form 
factors were set to zero. To take into account the influence of radiative corrections on the 
acceptance, the Photos package was used ([12]). The Dalitz plot density was then modified 
by the radiative corrections as calculated by Ginsberg [13]. This modification was made by 
randomly rejecting events according to the correct shape of the Dalitz plot. From 45 million 
generated MC events, 7.7 million were selected by the criteria in Sect. 3.2, corresponding to 
an acceptance of about 17%. 

Distributions for accepted Monte Carlo events weighted with a form factor slope of A + = 
0.029 are shown in figs. 2, 3 and 4 together with data, for the negative magnet polarity. 
As an example, the distributions of charged particle momenta in the Laboratory System 
(figs. 2, 3) and of the neutrino energy E* v in the centre-of-mass system, CMS, (fig. 4) are 
presented. The quantity E* is unambiguously defined. The upper part of the figures shows 
distributions of data, while the lower part shows the ratio of the experimental and MC 
spectra. Good agreement is seen for the Monte Carlo simulations with a form factor slope 
A , = 0.029. 
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events/1.4 GeV 

xlO 5 r 




100 



DATA/MC 



1.2 r- 
1.15 : 
1.1 : 




0.85 - I 

0.8 I" II' 

20 40 60 80 100 

Electron momentum (Lab System), GeV/c 

Fig. 2. Reconstructed electron momentum in the Laboratory System; data distribution and ratio 
data/MC 




Fig. 3. Reconstructed pion momentum in the Laboratory System; data distribution and ratio 
data/MC 

4 Data analysis 



4-1 Dalitz plot analysis 



For the Dalitz plot analysis, selected K e3 events were binned in a 19x19x19 three-dimensional 
array N(E*, q\/m\, q%/ml), where the indices 1 and 2 stand for the two possible solutions 
for the reconstructed kaon momentum and the bin sizes are 12 MeV for E* and 0.35 for the 
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events/2.3 MeV 

xlO 5 




Neutrino energy (CMS), GeV 



Fig. 4. Reconstructed neutrino energy E* in the Center of Mass System; data distribution and ratio 
data/MC 

q 2 /m\ bins. Two-dimensional arrays N(qf/m%, q%/m%) were used when pure V-A couplings 
were assumed. 

To derive the values of the form factors we used a maximum log-likelihood method [14] 
based on the following likelihood function which takes into account the finite size of the MC 
sample: 



InL = -2E(di In/, - /,) + In A, - A t )], (6) 

i i 

Here, di is the number of data events in the i'th bin, f\ is the predicted number of events, Oj is 
the actual number of MC events and is the expected number of MC events. The predicted 
number of events, fi, is a function of the three form factors parameters A+, fs, fr, and is 
obtained by weighting the MC expectation A{ for a flat Dalitz plot with the form factors 
according to Eqn.2. This procedure provides a weight hi for each bin. The MC expectation 
Ai for the bin % is obtained from the relation A^ = (di + aj)/(l + hi), which is valid at the 
maximum of the likelihood function [14], [15]. The predicted number f\ to be compared with 
the data di is then /j = A^ ■ hi. 

In order to test the goodness of fit, the variable 

X 2 gof = E (fi -di + di\nj)+Y^ (A -a l + ai In |-)] (7) 

was used [10], [16]. The minimization of the log-likelihood function, using MINUIT [17], gives 
the values of the fitted form factors. 
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iyf + (0)i 



Confidence levels: 

99 %,95 %,90 %,70 %,50 % 




0.03 



|f s /f + (0)l 

Fig. 5. Confidence level contours in the (|/s// + (0)|,|/ T /7+(0)|) plane at A+ = 0.0284. The out- 
ermost contour corresponds to 99% C.L., the innermost to 50% C.L. Only statistical errors are 
considered. 

DATA/MC 

1.6 i n 



1.5 - 




Fig. 6. Ratio of "diagonal" data and MC events (with A + = 0). The vertical axis corresponds to 
[/+(<? 2 )//+(0)] 2 - The curve corresponds to A+ = 0.0288 

4-2 Results 
Linear fits. 

We analysed the Dalitz plot under two hypotheses: either assuming pure V-A interaction 
or admitting additional scalar and tensor interactions. The results for the form factors ob- 
tained, assuming a linear q 2 dependence, are presented in Table 1. The statistical errors have 
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contributions from fluctuations of the data and of the MC events, as explained in sect. 4.1. 



ruiiii idLLui 


vdiue gluq ©GcLt. err. 


-v/ 2 /dot; 1 


A 


0284 + 0007 


8010/2915 


1 f<:/f^(0)\ 


015 + 2°°I 

U.U-L<J_Q 


3010/2915 


l/r//+(0)| 


o.o5±ig:S 


3010/2915 


A+ 


0.0288 ± 0.0004 


394/302 


A+ 


0.0280 ± 0.0019 


386/301 


a; 


0.0002 ± 0.0004 


386/301 



Table 1 

Results from the fits: 3-form factors case and 1-form factor case (last rows) 

The confidence level contours in the (|/s//+(0)|,|/t//+(0)|) plane for A + = 0.0284 are plotted 
in fig 5. The results are consistent with both fs and fx being zero at 10% C.L. taking only 
statistical errors into account. 

To visualize the main feature of the two-dimensional fit we display in fig. 6 the ratio of 
"diagonal" data and MC events (with A + = 0). "Diagonal" here means events for which both 
solutions for q 2 lie in one q 2 /m 2 -bin or they lie in neighbouring bins. The quantity plotted 
corresponds to the square of the q 2 variation of the vector form factor, [f + (q 2 ) / / + (0)] 2 = 
[1 + \ + q 2 /m 2 } 2 . The data are seen to be consistent with a linear increase of f+(q 2 ) with 
increasing q 2 . 

Our investigation of possible systematic errors showed that the biggest uncertainty for A + 
comes from the uncertainty in the kaon momentum spectrum. In order to determine the 
influence of this factor, the kaon momentum spectra derived from reconstructed K — > 7r + 7r~ 
and K — > 7r + 7r~7r° decays were implemented in the MC. The differences in the fitted form 
factor values using either of these spectra compared to the spectrum from diagonal K e3 
events was taken as a systematic error. This error amounts to 0.0007, and is 7 times larger 
than claimed in ref. [18]. We studied also the systematic errors induced by the various K e3 
selection cuts, by the momentum and energy calibrations, Dalitz plot distributions bin width, 
inefficiency of the muon veto system and the trigger, possible detector asymmetries and 
influence of accidental particles. 

The systematic uncertainties connected with our knowledge of the detector acceptance were 
evaluated by variations of the selection cuts in between values which changed the number of 
accepted events by up to 20%. The largest fluctuations in the form factor values were taken 
as systematic errors. 

Possible shifts in the measurement of p and E were investigated by comparing the value of 
the reconstructed kaon mass from K 37T and K 2n decays with the world average for the kaon 
mass and by comparing the central value of the E/p distribution of electrons with unity. 
These led to a rescaling of p and E by 1 • 10~ 3 and 4 • 10~ 3 respectively. The influence on the 
form factor values was studied by scaling systematically p and E by the quoted amounts. 
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Source 


A \ 
ZAA_|_ 


A 1 fs 1 

A l/ + (0)l 


A 1 /t 1 

A l/+(o)l 


only 


Kl spectrum 


±0.000o0 


i n am 

±0.0U1 


in nncr 
±0.005 


±0.000 (U 


Ljeom. accept. 


±0.00050 


±0.0U/ 


in ni r 
±0.015 


in r\r\r\Ar\ 
±0.00040 


Pmin 


±0.00025 


±0.0U4 


i n ni n 
±0.010 


in nnni c 
±0.00015 


L>lkr 


±U.UUU4b 


_i_n nn/i 
±U.UU4 


in nnc^ 
±U.UUO 


in nnnoc^ 
±U.UUU2o 


h/p 


i n nnnor; 

±0.00035 


±0.0U2 


in n-i n 
±0.010 


in nnnoc; 
±0.00035 




±0.00020 


±0.0Uo 


in nncr 
±0.005 


in nnni n 
±0.00010 


E,p scaling 


±U.UUU2U 


±U.UU1 


i n nnr 

±U.UUO 


i n finnon 
±U.UUU2U 


MU V men. 


±0.00020 


±0.002 


i n nnc; 

±0.005 


i n nnnin 

±0.00020 


Trigger ineff. 


_i_n nnm k 

±U.UUU1D 


_i_n nno 
±u.uuz 


_i_n nn^ 
±u.uuo 


_i_n nnno^ 


accidentals 


±0.00030 


±0.001 


±0.005 


±0.00025 


bin width 


±0.00040 


±0.005 


±0.010 


±0.00010 


TOTAL 


±0.0013 


±0.012 


±0.03 


±0.0011 



Table 2 

Systematic uncertainties: 3-form factors case and 1-form factor case (right column) 

The inefficiency of the MUV during this run was measured to be below 2 %. In order to 
investigate the influence of the MUV inefficiency we introduced an extra inefficiency of 2 % 
in the experimental data and compared the results. The difference was taken as a systematic 
error. 

The trigger inefficiency was evaluated by using the control trigger and was found to be (1.9 ± 
0.1)%. The influence of this source of systematic uncertainty on the results was estimated 
as follows. The results obtained with the main trigger data for which the control trigger was 
set as well were compared to the results obtained with the control trigger data which are 
strongly correlated. The differences are due to the 2% events escaping the main trigger, and 
indicate the uncertainty due to this inefficiency. 

To estimate the possible influence of accidental particles we used the Overlay MC technique, 
applied earlier in the e'/e analysis [11], in which randomly recorded experimental signals with 
rate proportional to the beam intensity were mixed with each MC event. The uncertainty 
due to this effect was estimated by comparing results obtained with the standard MC with 
results obtained with the Overlay MC. 

We also varied the number of bins (from 10 to 25 in each axis) in the distributions 
N(E v ,qf,q%) and N{q\,q%). The change in the values of the fit parameters was considered 
as a systematic error. 

The individual systematic uncertainties, and the total systematic error obtained by combin- 
ing the individual errors in quadrature, are summarized in table 2. 

As a cross check, correlations between magnetic field polarities, positive and negative tracks, 
particles and anti-particles, geometrical positions of the tracks were explored. We did not 
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find significant changes in results coming from these sources. The effect of the uncertainties 
of the radiative corrections [13] was found to have a negligible effect on the results. 

Quadratic and dipole fits. 

As an alternative to these linear slope fits, we also investigated two other forms of the 
q 2 dependence: a quadratic form, admitting a finite A' + as a free parameter, and a pole- 
dominance form factor of the form given in Sect.l, with the mass of the vector meson pole 
mass My as the only free parameter. The quadratic fit gave the values of the two slope 
coefficients, A+ = (28.0 ± 1.9 ± 1.5) • 1(T 3 and A' + = (0.2 ± 0.4 ± 0.2) • 10~ 3 , consistent with 
no quadratic term but also consistent with a dipole form factor where the quadratic term 
of the Taylor expansion corresponds to X' + = \+ = 0.7 ■ 10~ 3 . The dipole fit with the vector 
meson mass as free parameter gave My = 859 ± 18 MeV, fully consistent with the mass of 
the AT* (892) meson. 



5 Conclusions 



Using 5.6 million K e3 events recorded by the NA48 detector, we tested the validity of the 
V-A weak interaction and measured the q 2 dependence of the vector form factor with high 
precision. Admitting scalar and tensor couplings in addition to V-A couplings , we obtain 
the following form factor values : 

|-^y| = o.oi5±8:8S5 ±0.012 

1^1 = 0.0518$ ±0.03 



and the linear slope in the vector form factor A + = (28.4 ± 0.7 ± 1.3) • 10 3 . The 90% C.L. 
upper limits on the scalar and tensor form factors are I I < 0.041 and I /J n . 1 < 0.12. We 
therefore see no evidence for scalar or tensor couplings, not supporting the results of ref . [6] . 



Assuming only V-A currents and a linear q 2 dependence the A + value is A + = (28.8±0.4±1.1)- 
10~ 3 . If in addition we admit quadratic terms , the results are A + = (28.0±1.9±1.5)-10~ 3 and 
A' + = (0.2 ± 0.4 ± 0.2) • 10~ 3 . We therefore find no evidence for a quadratic term A' + different 
from zero. However, the size of a quadratic term is consistent with a Taylor expansion of a 
pole-dominance form factor. Such a dipole form factor is in good agreement with the data, 
with a pole mass of 859 ± 18 MeV. 

Fig. 7 shows the linear slope. 

Our results agree well with the V-A coupling of the weak interaction and with a pole form of 
the form factor, in agreement with the linear approximations tested in former experiments 
[9], [18], [19], [20], [21], [22], [23]. Our result for the quadratic term in the vector form factor 
is at variance with the result of ref. [18] though it agrees with the pole model form factor. 
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Fig. 7. Experimental results for the slope A + of the vector form factor. The bold error bars are the 
statistical errors, while the thin line indicates the total uncertainty including the systematic error. 
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